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Abstract 

The  objective  of  this  work  is  to  examine  the  reactant  gas  transport  and  the  cell  performance  of  a  proton  exchange  membrane  fuel  cell  (PEMFC) 
with  a  tapered  flow  channel  design.  It  is  expected  that,  with  the  reduction  in  the  channel  depth  along  the  streamwise  direction,  the  reactant  fuel  gas 
in  the  tapered  channel  can  be  accelerated  as  well  as  forced  into  the  gas  diffuser  layer  to  enhance  the  electrochemical  reaction  and  thus  augment 
the  cell  performance.  The  effects  of  liquid  water  formation  on  the  reactant  gas  transport  are  taken  into  account  in  the  present  study.  Numerical 
predictions  show  that  the  cell  performance  can  be  enhanced  with  the  fuel  channel  tapered,  and  the  enhancement  is  more  noticeable  at  a  lower 
voltage.  The  results  also  reveal  that  the  liquid  water  effect  in  general  influences  the  cell  performance  and  the  effect  becomes  significant  at  lower 
voltages. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  is  a  class  of  electro-chemical  devices  that  realizes 
the  direct  conversion  of  chemical  energy  of  reactants  to  electric 
energy.  With  high  efficiency  and  high  environment  compatibil¬ 
ity,  the  proton  exchange  membrane  (PEM)  fuel  cell  has  become  a 
potential  candidate  of  the  power  source  in  the  near  future.  There 
have  appeared  a  number  of  studies  on  the  transport  phenomena 
and  the  cell  performance  of  the  PEM  fuel  cells  in  the  open  liter¬ 
ature.  Springer  et  al.  [1]  proposed  a  one-dimensional  isothermal 
model  for  the  PEM  fuel  cells.  In  their  model,  water  diffu¬ 
sion  coefficient,  electro-osmotic  drag  coefficient,  water  sorption 
isotherms  and  the  membrane  conductivities  were  all  assumed  to 
be  functions  of  membrane  water  content.  It  is  found  that  the 
membrane  resistance  increases  with  the  operating  current  den¬ 
sity  and,  at  typical  operating  conditions,  the  net  water  flux  ratio 
is  much  less  than  that  with  a  fully  hydrated  membrane.  It  is  also 
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disclosed  that  the  resistance  can  be  reduced  for  the  membrane 
thin  enough. 

Nguyen  and  White  [2]  proposed  a  model  with  the  effects 
of  electro-osmosis,  diffusion  of  water,  heat  transfer  from  solid 
to  gas  and  latent  heat  of  water  evaporation/condensation  to 
study  the  water  and  thermal  management  of  PEM  fuel  cells. 
They  found  that  the  ohmic  loss  is  considerable  at  high  current 
density  and  the  voltage  loss  is  twice  amount  of  that  on  the  cath¬ 
ode  electrode.  The  back  diffusion  of  the  water  from  cathode  to 
anode  is  not  enough  to  keep  the  membrane  hydrated,  therefore, 
the  reactant  gas  at  the  anode  needs  to  be  humidified.  Yi  and 
Nguyen  [3]  developed  an  along-the-channel  model  for  evalu¬ 
ating  the  effects  of  various  designs  and  operating  parameters 
on  the  performance  of  a  PEM  fuel  cell.  The  results  showed 
that  the  humidification  of  the  anode  gas  is  required  to  improve 
the  conductivity  of  the  membrane,  and  the  liquid  injection  and 
higher  humidification  temperature  can  enhance  the  cell  perfor¬ 
mance.  Additionally,  applying  higher  cathode  gas  pressure  is 
helpful  to  replenish  the  water  loss  due  to  electro-osmosis  and, 
thereby,  results  in  more  conductive  membrane  and  better  cell 
performance. 
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Nomenclature 

C  mass  fraction 

Cp  quadratic  drag  factor 

D  mass  diffusivity  (m2  s-1) 

Du+  proton  diffusivity  in  membrane 

(4.5  x  10-9  atom2  s-1)  (Gurau  et  al.  [19]) 

Z)r2  hydrogen  diffusivity  (2.63  x  10-6  m2  s-1) 

(Um  and  Wang  [9]) 

Do2  oxygen  diffusivity  (5.2197  x  10-6m2s-1) 

(Um  and  Wang  [9]) 

F  Faraday  constant  (96487  C  mol- 1 ) 

I  current  density  (A  m-2) 

jQCf  reference  exchange  current  density  at  anode  side 

(5  x  102  Am-3)  (Gurau  et  al.  [19]) 
kc  condensation  rate  constant  (s-1)  (Mazumder  and 
Cole  [6]) 

kQ  evaporation  rate  constant  (atm- 1  s  - 1 )  (Mazumder 

and  Cole  [6]) 

kv  permeability  (m2) 

N  number  of  baffles 

P  pressure  (Pa) 

R  universal  gas  constant  (8.314  J  mol-1  K-1) 

Rch  ratio  of  the  outlet  channel  depth  to  inlet  channel 
depth 

.v  saturated  rate 

S  source  term 

T  temperature  (K) 

u  velocity  in  the  v-direction  (ms-1) 

v  velocity  in  the y-direction  (ms-1) 

V  operating  voltage  (V) 

w  oxygen  mass  fraction 

vw  mole  fraction  of  water 

Zf  spices  valence 

Greek  letters 

a  transfer  coefficient  for  the  reaction  (0.5  on  anode 

side,  0.5  on  cathode  side) 
s  porosity 

0  surface  overpotential  (V) 

p  density  (kg  cm-  3 ) 

a  ionic  conductivity  of  the  ionomer  (£2-1  m-1) 

r  tortuosity  of  diffusion  layer 

v  viscosity  of  flow  (kg  m-1  s-1) 

uair  air  viscosity  (1.58  x  10-5  kgm-1  s-1) 

(Gurau  et  al.  [19]) 

URoO  water  viscosity  (8.91  x  10  4 kgm  1  s  1 ) 

(Gurau  et  al.  [19]) 

0  ionomer  phase  potential  (V) 

0  phase  potential  (V) 

Subscripts 
a  anode 

c  cathode 

d  diffuser  layer 

eff  effective 


k  kth  composition  of  fuel  reactant 

L  liquid  water 

m  membrane 

sat  saturated 

u  v-direct 

v  y-direct 

Superscripts 
ref  reference 


Ge  and  Yi  [4]  developed  a  two-dimension  model  to  investi¬ 
gate  the  effects  of  operation  conditions  and  membrane  thickness 
on  the  water  transport.  In  their  study,  the  change  in  effec¬ 
tive  porosity  with  liquid  water  was  considered  to  simplify  the 
model  of  the  two-phase  transport  phenomena  in  porous  layers. 
The  results  revealed  that  the  cell  performance  can  be  enhanced 
by  increasing  the  cell  temperature.  Mazumder  and  Cole  [5,6] 
numerically  investigated  the  cell  performance  of  PEM  fuel  cells 
with  consideration  of  the  liquid  water  effects.  The  results  showed 
that  as  the  current  density  gets  up  to  a  certain  critical  value,  the 
percentage  of  the  porous  media  filled  with  liquid  water  may 
exceed  50%,  especially  at  the  cathode  side. 

The  flow  field  designs  in  bipolar  plates  are  important  in  the 
development  of  the  PEM  fuel  cells.  Nguyen  [7]  proposed  a  new 
flow  field  design  to  enhance  the  liquid  water  transport  out  of 
the  gas  diffuser  layer  (GDL).  This  design  works  by  converting 
the  transport  of  reactant/product  gases  to/from  the  catalyst  layers 
from  diffusion  mechanism  to  a  convection  mechanism.  Um  et  al. 
[8],  and  Um  and  Wang  [9]  developed  a  multi-dimensional  model 
to  study  the  electrochemical  kinetics,  current  distributions,  fuel 
and  oxidant  flows,  and  multi-component  transports  in  a  PEMFC 
with  the  interdigitated  flow  field.  To  understand  the  effects  of 
the  interdigitated  gas  distributor  on  the  performance  of  these  fuel 
cell  electrodes,  Yi  and  Nguyen  [10]  employed  a  two-dimensional 
isothermal  model  of  a  porous  electrode  to  simulate  the  hydrody¬ 
namics  of  gas  flow  through  the  pore  volume  of  the  electrode  of 
a  PEMFC.  It  is  concluded  that,  from  the  predictions  of  the  PEM 
fuel  cells  with  interdigitated  flow  field,  the  higher  gas  flow  rate 
through  the  electrode  improves  the  electrode  performance. 

He  et  al.  [11]  studied  the  influences  of  electrode  and  flow 
field  design  on  the  performance  of  a  PEM  fuel  cell  with  a 
half-cell  model.  It  was  found  that  higher  differential  pressure 
between  inlet  and  outlet  channels  would  enhance  the  electrode 
performance.  Increasing  the  electrode  thickness  is  equivalent 
to  increasing  the  diameter  of  a  pipe  in  a  fluid  flow  system 
and  the  gas  bypass  effect  becomes  more  significant.  Dohle 
et  al.  [12]  examined  the  effects  of  flow  fields  on  the  power 
density  and  efficiency  of  fuel  cell  systems  by  conducting 
experiments  and  numerical  simulations.  In  their  work,  the 
see-through  bipolar  plate  and  additional  resistance  matrix  (with 
simulation)  are  employed  to  measure  the  flow  propagation 
velocity  and  flow  distribution  with  different  permeability. 
The  electrochemical  behaviors  of  a  direct  methanol  fuel  cell 
(DMFC)  with  serpentine  or  interdigitated  flow  field  at  both 


80 


H.C.  Liu  et  al.  /  Journal  of  Power  Sources  158  (2006)  78-87 


cathode  and  anode  were  investigated  by  Arico  et  al.  [13]. 
The  experimental  results  showed  that  the  interdigitated  design 
significantly  enhances  the  mass -transport  inside  a  DMFC  and 
achieves  the  higher  maximum  power  outputs  compared  to  the 
classical  serpentine  geometry.  The  larger  methanol  permeation 
through  the  electrolyte  determines  both  lower  voltage  and  fuel 
efficiencies  in  the  activation-controlled  region  with  respect  to 
the  serpentine  design.  Recently,  Yan  et  al.  [14,15]  studied  the 
effects  of  flow  channel  designs  on  the  steady  or  dynamic  gas 
transport  and  cell  performance  of  PEM  fuel  cells.  They  found 
that  the  flow  distributor  geometry  has  a  significant  influence 
on  the  cell  performance.  In  addition,  they  indicated  that  for  a 
larger  channel  depth  fraction,  a  faster  dynamic  response  as  well 
as  better  cell  performance  can  be  attained. 

Recently,  Yan  and  co-workers  [16,17]  proposed  a  new  flow 
field  with  baffles  in  the  tandem  array  to  augment  the  reactant 
transport  and  cell  performance  of  a  proton  exchange  membrane 
fuel  cell.  They  found  that  the  reactant  transport  and  cell  perfor¬ 
mance  both  can  be  enhanced  by  the  presence  of  the  baffles  in 
the  flow  channel  of  the  bipolar  plate,  especially  at  the  operating 
conditions  of  low  voltage.  The  beneficial  baffle  effects  become 
increasingly  remarkable  with  increasing  width  and/or  number 
of  baffles  in  the  tandem  array.  Soong  et  al.  [18]  proposed  a  rel¬ 
atively  novel  configuration  of  partially  blocked  fuel  channels. 
The  blockage  effects  with  various  gap  ratios  and  baffle  num¬ 


Y 


X 


bers  on  the  reactant  gas  transport  and  the  pressure  loss  across 
the  channel  length  were  emphasized.  It  was  found  that  reducing 
gap  size  and/or  increasing  baffle  number  to  enhance  the  reactant 
gas  transport  may  result  in  a  penalty  of  high  pressure-loss.  Very 
high  pressure-loss  can  be  generated  due  to  high  flow  resistance 
at  a  low  gap  ratio  combined  with  a  low  GDL  porosity.  They 
also  proposed  reasonable  combinations  of  the  influential  design 
parameters. 

From  the  literature  survey  presented  above,  it  is  found  that  the 
forced  convection  of  fuel  due  to  channel  blockage  always  leads 
to  a  pressure  loss.  To  improve  the  performance  of  PEM  fuel  cells 
but  with  the  penalty  of  high  pressure-loss  alleviated  is  one  of  the 
significant  issues  in  development  of  the  fuel  cell  technology.  The 
partially  blocked  channel  design  studied  previously  is  one  of  the 
trials  and  other  possibilities  are  still  open.  The  objective  of  the 
present  study  is  to  investigate  a  new  design  of  fuel  channel  with 
tapered  configuration  in  main  flow  direction.  A  two-dimensional 
numerical  model  of  PEMFC  with  tapered  flow  channel  design 
at  cathode  is  developed  to  employee  to  perform  the  exploration 
of  the  resultant  transport  phenomena  and  cell  performance. 

2.  Analysis 

The  schematic  diagram  of  the  physical  system  under  consid¬ 
eration  is  shown  in  Fig.  1,  in  which  the  PEMFC  consists  of  a 


w 

vv  out  N — 


Fig.  1.  Schematic  diagram  of  2D  PEMFC  model. 
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Table  1 

Detailed  expressions  of  the  source  terms  in  the  governing  equations 
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proton  exchange  membrane  in  the  middle,  and  two  catalyst  lay¬ 
ers,  two  gas  diffuser  layers,  and  two  flow  channels  each  at  the 
anode/cathode  side.  It  is  seen  that  the  flow  channel  depth  at  the 
cathode  side  reduces  gradually  along  the  main  flow  direction.  At 
the  cathode  side,  the  reduction  in  flow  area  of  the  fuel  channel 
will  force  the  reactant  gas  transfer  to  GDL  and  catalyst  layer 
as  well  as  increase  the  forced  convection  in  the  flow  direction. 
Firstly,  the  forced  convection  in  diffuser  layer  would  generate 
the  shear  stress  to  remove  the  liquid  water  and  heat  to  reduce 
the  problem  of  water  and  thermal  management  in  the  PEM  fuel 
cell.  Secondly,  the  reactant  gas  is  easier  to  react  in  the  catalyst 
layer,  which  in  turn,  causes  a  better  cell  performance.  There¬ 
fore,  it  is  appealing  to  have  a  numerical  study  on  the  effects  of 
tapered  flow  channel  design  on  the  cell  performance  of  PEM  fuel 
cells.  To  simplify  the  analysis  in  the  present  study,  the  following 
assumptions  are  invoked:  (1)  the  gas  mixtures  are  considered  to 
be  perfect  gases;  (2)  the  gas  flow  is  steady  and  laminar  every¬ 
where  in  the  cell;  (3)  the  gas  diffuser  layers,  catalyst  layers,  and 
the  membrane  are  considered  to  be  isotropic  porous  media;  (4) 
the  volume  of  liquid  water  in  the  PEM  fuel  cell  is  negligible; 
and  (5)  the  heat  generated  by  chemical  reaction  is  negligible. 

With  the  above  assumptions,  the  gas  transport  equations  for 
the  two-dimensional  PEM  fuel  cell  can  be  described  as  follows. 

Continuity  equation: 


du  dv 
dx  +  dy 


Momentum  equations: 
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The  source  terms,  Su,  Sv  and  Sc,  in  the  momentum  and  species 
concentration  equations  are  listed  in  Table  1 .  Among  them,  the 
source  terms  Su  and  Sv  account  for  the  Darcy’s  drag  forces  in  the 
v-  and  y-directions,  respectively.  The  parameter  £eff  stands  for 
the  effective  porosity,  Cp  the  quadratic  drag  factor,  kv  the  per¬ 
meability,  and  D^eff  the  effective  diffusion  coefficient  of  the  kth 
composition  of  fuel  reactant  determined  by  Bruggerman  model 
from  He  et  al.  [11].  In  Eq.  (4),  Sc  and  Sl  are  the  source  term 
of  chemical  reaction  and  liquid  water  in  the  species  concentra¬ 
tion  equation,  and  Sc  is  different  for  the  reactant  gases,  e.g.  Sc 
is  -7a/2FCtotai,a  f°r  hydrogen,  -  jc/4FCtotal,c  for  oxygen,  and 
— 7c/2FCtotai,c  f°r  water  vapor.  The  parameters  ja  and  jc  indicate 
the  current  densities  at  the  anode  and  cathode  sides,  respectively, 
and  can  be  described  by  the  following  Butler- Volmer  equations 
[19]: 


Q(aaF/RT)r] 
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e(ac  F/RT)r] 


Q(aaF/RT)r] 


1 

e(ac  F/RT)r] 


In  the  above  equations,  A  j gef  is  the  reference  exchange  current 
density,  and  ac  are  the  electric  charge  transport  rates  in  anode 
and  cathode  catalyst  layers,  rj  is  the  overpotential,  F  the  Fara¬ 
day’s  constant,  R  the  gas  constant  and  T  the  temperature  of  the 
fuel  cell.  The  parameters  used  are  listed  in  Table  2. 

In  order  to  evaluate  the  distributions  of  the  local  current  den¬ 
sity,  the  following  phase  potential  equation  is  employed, 
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As  for  the  transport  equation  of  the  proton  in  the  membrane, 
one  has 
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Table  2 

Parameters  used  in  this  work  (Gurau  et  al.  [19]) 


Porosity  of  diffuser  layer 
Permeability  of  diffuser  layer 
Porosity  of  catalyst  layer 
Permeability  of  catalyst  layer 
Porosity  of  membrane 
Permeability  of  membrane 
Baffle  width 
Channel  length 
Channel  depth  (anode  side) 
Diffuser  layer  thickness 
Catalyst  layer  thickness 
Membrane  thickness 
Operation  temperature 
Operation  pressure 
Inlet  average  velocity  (i>o) 
Relative  humidity(both  anode  and 
Cell  temperature 
Anode  fuel 
Cathode  fuel 


0.4 

1.76  x  10-11  m2 
0.4 

1.76  x  10-11  m2 
0.28 

1.58  x  10-18  m2 
0.000762  m 
0.0762  m 
0.000762  m 
0.000254  m 
0.0000287  m 
0.00023  m 
353  K 
1  atm 
1.76ms_1 

cathode)  100% 

80°C 

H2 

Air 


where  0  is  the  phase  potential  and  crm  is  the  electric  conductiv¬ 
ity  of  the  membrane  which  can  be  calculated  by  the  formulas 
developed  by  Springer  et  al.  [1]: 


crm(T)  =  o+ex  P 


where  the  reference  electric  conductivity  is 

oief  =  0.005 139A-  0.00326  (10) 


0.043  +  17.81a  —  39. 85a2 +  36.0a3  0  <  a  <  1 
14  +  1.4(a  —  1)  1  <  a  <  3 


(11) 


In  Eq.  (1 1),  a  =  M^o^V^sat  Is  the  activity  of  water  and  the  satu¬ 
ration  at  the  pressure  of  water  Psat  is  calculated  by  the  following 
equation  [1], 

p  t  _  10—2.1794+0.029537— 9. 1837x  10-5r2  +  1.4454x  10-7  q2) 


Using  the  following  relations  between  the  phase  potential  0 
and  current  density  i, 
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Eq.  (7)  can  then  be  rewritten  as 
dix  di 
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-y 


dx  dy 


=  S 


(15) 


The  water  management  is  very  important  to  the  performance 
of  PEM  fuel  cells.  In  practice,  the  liquid  water  in  the  porous 
material  could  occupy  the  pore  and  thus  reduce  the  ability  of  the 
gas  transport  in  the  porous  layers.  For  simplicity,  only  the  effects 
of  liquid  water  on  the  effective  porosity  and  mass  diffusivity  in 
the  gas  diffuser  and  catalyst  layers  are  taken  into  account.  As 
the  partial  pressure  of  water  exceeds  the  saturation  pressure  of 
water,  the  following  equation  is  used  to  determine  the  sources 


term  of  liquid  water,  S l  [6], 


P«,_o  -  ft.),  if  Pn,o  >  />«, 

ke£efU(PSat  -  ^H2C),  if  ^H20  <  ^sat 

(16) 


where  the  kc  and  kQ  are,  respectively,  the  constants  of  water 
condensation  and  evaporation  [6],  and  the  PSat  is  evaluated  by 
using  Eq.  (12). 

The  saturation  rate  s  is  employed  to  express  the  ratio  of  the 
volume  of  pore  occupied  by  the  liquid  water  to  the  volume  of 
pore  in  the  porous  medium,  viz. 

the  volume  of  pore  occupied  by  liquid  in  porous  medium 
v  =  - 

total  pore  volume  in  porous  medium 

(17) 


The  effective  porosity  of  porous  medium  is  modified  to 
account  for  the  liquid  water  effect, 


£eff  =  e(l  -  s) 


(18) 


Boundary  conditions  for  the  solutions  of  the  transport  equa¬ 
tions  at  the  interfaces  between  different  layers  of  the  same 
domain  are  not  required.  At  the  gas  channel  outlet,  the  fully 
developed  flow  conditions  are  assumed, 

dv  dCk 


u  =  —  =  -  =  0 

dy  dy 


(19) 


The  boundary  conditions  at  the  channel  walls  and  the  plate 
collectors  are 


dCk 

u  =  v  =  - =  0 

dx 


(20) 


At  the  interfaces  between  the  gas  diffusers  and  the  gas  chan¬ 
nels,  the  following  boundary  conditions  are  used, 


£&ff,X+ 


£&ff,X+ 


dv 
dx 

dCk 

dx 


dv 


x=X+ 


dx 


vx=X+  —  Vx=X- 


(21) 


x=X~ 

dCk 


x=X 


+ 


dx 


Ck,x=X+  ~  Ck,x=X~  (22) 


1 


The  similar  conditions  are  employed  for  the  interface  between 
the  gas  diffuser  and  the  catalyst  layer  can  be  expressed  as  fol¬ 
lows, 


£Qff,X+ 


dv 

dx 


dv 

=  £eflX~  ^ 


vx=x+  =  vx=x~  (23) 


x=X- 


£Qff,X+ 


dCt 

dx 


dCk 

~  £eff,X-  -X- 
x=X+  °x 


x=X- 


Ck,x=x+  —  Ck,x=x~  ■  (24) 

The  interface  between  the  catalyst  layer  and  membrane  is 


dv 

dx 


=  0, 


CL 


9Q 

dx 


=  0. 


(25) 


CL 


The  boundary  conditions  for  the  phase  potential  at  the  inter¬ 
face  between  the  catalyst  layer  and  the  membrane  are  0  =  0 
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on  the  anode  side  and  d@/dX=0  on  the  cathode  side.  Since  the 
phase  potential  is  a  linear  distribution  in  the  membrane,  the  phase 
potential  boundary  condition  can  be  written  as  d@/dY=  0. 

3.  Numerical  method 

The  solution  to  the  governing  equations  is  performed  by 
employing  a  finite  volume  scheme  with  the  model  domain 
divided  into  a  number  of  cells  as  control  volumes.  The  gov¬ 
erning  equations  are  numerically  integrated  over  each  of  these 
computational  cells  or  control  volumes.  The  method  exploits 
a  collocated  cell-centered  variable  arrangement  with  the  local 
or  cell-averaged  values  of  the  physical  quantities  evaluated  and 
stored  at  each  cell  center. 

The  governing  equations  can  be  expressed  in  the  form  of  a 
generalized  convection-diffusion  type  of  transport  equation: 

V  (/ ouct>  -  r^V<j>)  =  S<p  (26) 

where  0  denotes  the  general  dependent  variable,  r 0  the 
exchange  coefficient,  the  source  term,  u  the  velocity  vector, 
and  p  the  density.  With  the  discretization  of  the  governing  equa¬ 
tions,  the  coupled  finite-difference  equations  can  be  expressed 
in  the  form  of 

apcpp  =  <2e4>e  +  aw^w  +  un&n  +  as&s  +  Sf  (27) 

where  <fip  is  the  value  of  0  at  the  current  point  P,  0#  •  •  •  0s  stand 
for  the  values  of  the  grid  points  adjacent  to  the  point  P,  and 
ap  •  •  •  as  are  known  as  the  link  coefficients. 

In  this  work,  the  non-uniform  grid  system  of  72  x  51  points 
is  used.  To  examine  the  grid  independence  of  the  predictions, 
three  grid  systems  shown  in  Fig.  2  are  considered  and  their  influ¬ 
ences  on  the  prediction  of  local  current  density  distribution  for 
a  typical  case  are  presented  in  Fig.  3.  It  is  found  that  the  maxi¬ 
mum  deviation  among  the  computations  on  the  grids  of  37  x  26, 
72  x  5 1  and  1 13  x  73  are  less  than  3%  and  the  results  on  72  x  51 


Fig.  3.  Comparison  of  predictions  on  the  three  different  grid  systems. 

and  113  x  73  grids  are  quite  close.  Therefore,  the  grid  system  of 
72  x  51  points  seems  sufficient  to  resolve  the  behaviors  of  the 
reactant  gas  transport  in  the  present  PEM  fuel  cell  model. 

4.  Results  and  discussion 

To  examine  the  effects  of  the  tapered  flow  channel,  a  param¬ 
eter,  7?ch,  characterizing  the  reduction  in  flow  channel  depth  is 
defined  as  the  ratio  of  channel  depth  at  the  outlet  to  that  at  the 
inlet.  In  this  work,  the  channel  depth  at  the  inlet  is  0.762  mm. 
Therefore,  the  channel  depths  at  exit  are,  respectively,  0.762  mm 
for7?ch  =  1,0.381  mmfor7?ch  =  0.5  and  0.0762  mm  for  Rch  =  0.1. 
The  effects  of  Rc h  on  the  I-V  and  I-P  curves  are  presented  in 
Fig.  4.  An  overall  inspection  on  Fig.  4  discloses  that,  under  the 
conditions  of  the  high  operating  voltages  or  low  current  densi¬ 
ties,  the  7?ch  effects  on  the  cell  performance  are  negligibly  small. 
At  the  low  operating  voltages,  however,  the  variation  in  Rc h  sig¬ 
nificantly  alters  the  performance  curves.  Relative  to  the  PEMFC 
with  the  channel  of  uniform  width  (Rc h  =  1),  the  improvement 


O  Fuel  outlet 


Fuel  outlet 


Cathode 


Fuel  inlet 


Fuel  inlet 


Fig.  2.  Mesh  structure  of  three  different  grid  systems  of  grid  points:  (a)  113  x  73;  (b)  72  x  51;  and  (c)  37  x  26. 
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Fig.  4.  Effects  of  Rc h  on  the  cell  performance:  (a)  I-V curves;  and  (b)  I-P  curves 
(note:  P  stands  for  power  in  this  plot). 


0  0.01  0.02  0.03  0.04  0.05  0.06  0.07 


(b)  Y(m) 

Fig.  5.  Effects  of  Rc h  on  the  local  current  density  distributions  at  (a)  V=  0.7  V; 
and  (b)  V=0.2V. 


in  the  cell  performance  with  the  tapered  flow  channel  is  demon¬ 
strated.  It  is  clearly  observed  that  the  smaller  the  Rc h,  the  higher 
the  limiting  current  density  is.  The  liquid  water  effects  are  also 
shown  in  Fig.  4.  At  high  operating  voltages,  the  deviations  in 
the  cell  performance  between  the  results  with  and  without  con¬ 
sideration  of  liquid  water  effects  are  negligible.  This  means  that 
the  fuel  transport  in  the  PEM  fuel  cell  can  be  safely  modeled 
as  a  single-phase  flow.  While  at  low  voltages,  the  liquid  water 
effects  on  the  cell  performance  become  remarkable  and  cannot 
be  neglected  in  the  modeling.  This  confirms  the  fact  that  the  mass 
transport  at  lower  voltages  is  more  significant  and,  in  turn,  more 
water  is  generated  in  the  catalyst  layer  of  the  cathode  side.  The 
two-phase  flow  effects  should  be  considered  at  low  voltages. 

The  7?ch  effects  on  the  local  current  density  distributions  at 
V=0.7  and  0.2  V  are,  respectively,  presented  in  Fig.  5(a)  and 
(b),  in  which  the  predictions  with  and  without  consideration 
of  liquid  water  are  compared  with  the  corresponding  cases  of 
the  results  without  liquid  water  effects.  By  comparing  the  ordi¬ 
nates  in  Fig.  5(a)  and  (b),  it  indicates  that  the  variations  in  the 
local  current  density  along  the  flow  direction  are  relatively  more 
noticeable  in  the  cases  of  a  low  operating  voltage  (V=0.2V). 
The  electrochemical  reaction  in  the  catalyst  layer  is  expected 
to  be  stronger  at  a  lower  operating  voltage  or  a  higher  operat¬ 
ing  current  density.  Therefore,  the  reactant  oxygen  consumes  at 
a  higher  rate  as  it  flows  along  the  main  flow  direction,  which 


Fig.  6.  Effects  of  Rc h  on  the  local  mass  fraction  distributions  along  the  interface 
between  the  GDL  and  catalyst  layer  at  V=  0.7  V:  (a)  oxygen;  and  (b)  hydrogen. 
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in  turn  causes  the  substantial  variations  in  the  local  current  den¬ 
sity.  Fig.  5  also  discloses  that  the  better  cell  performance  is  noted 
for  the  fuel  channel  with  a  smaller  Rc h.  In  addition,  the  results 
without  consideration  of  liquid  water  over-predict  the  cell  per¬ 
formance. 

In  the  study  of  reactant  gas  transport  in  PEM  fuel  cells,  under¬ 
standing  of  the  detailed  distributions  of  the  fuel  gases  at  the 
interface  between  the  gas  diffuser  layer  and  catalyst  layer  are 
important  to  the  design  of  a  PEM  fuel  cell.  To  this  end,  the 
Rch  effects  on  the  local  distributions  of  oxygen  and  hydrogen 
mass  fractions  at  the  GDL-catalyst  interface  at  operating  volt¬ 
age  V=0.2  V  are  presented  in  Fig.  6.  The  oxygen  mass  fraction 
on  the  cathode  side  decreases  along  the  axial  location  Y  due  to 
the  chemical  reaction  of  the  oxygen.  Results  in  Fig.  6(a)  show 
that,  near  the  entrance  (T <  0.01  m),  the  oxygen  fraction  is  lower 
for  the  channel  with  a  lower  Rc h  (=0.1),  in  which  more  reac¬ 
tant  gas  in  the  downstream  part  ( Y >  0.01  m)  of  the  channel  can 
be  forced  into  the  GDL  and  catalyst  layers  and  thus  result  in  a 
higher  oxygen  fraction.  It  is  also  found  in  Fig.  6(a)  that  the  oxy¬ 
gen  faction  is  relatively  lower  for  the  case  with  consideration 
of  the  liquid  water  effects.  This  is  due  to  the  reduction  in  the 
oxygen  diffusion  with  the  appearance  of  the  liquid  water,  which 
in  turn  causes  a  lower  oxygen  fraction.  As  for  the  Rc h  effects  on 
the  anode  reactant  gas,  as  shown  in  Fig.  6(b),  the  hydrogen  mass 
faction  distributions  are  only  slightly  affected  by  the  variations 


Fig.  7.  Effects  of  Rc h  on  the  local  liquid  water  distributions  along  the  interface 
between  the  GDL  and  catalyst  layer  at:  (a)  V  =  0.7  V;  and  (b)  V  =  0.2  V. 


in  Rc h.  This  confirms  the  general  concept  of  the  negligible  effect 
of  the  cathode  flow  field  on  the  anode- side  transport  phenomena. 

Thermal  and  water  management  is  another  important  issue 
pertaining  to  the  performance  of  PEM  fuel  cells.  To  explore  the 
effect  of  the  tapered  flow  channel  on  the  water  content  distri¬ 
butions,  the  liquid  water  saturations  at  the  GDL-catalyst  layer 
interface  for  different  Rc h  are  predicted  and  shown  in  Fig.  7,  in 
which  the  solid  and  dash  lines  denote  the  results  on  the  cath¬ 
ode  and  anode  sides,  respectively.  With  the  same  Rc h,  the  liquid 
water  saturation  S  is  higher  in  a  case  of  lower  operating  voltage. 
This  can  be  made  plausible  by  noting  that,  for  a  lower  operating 
voltage  (or  a  higher  operating  current  density),  the  electrochem¬ 
ical  reaction  is  faster  and  thus  generates  more  liquid  water.  The 
noticeable  Rc h  effects  on  the  liquid  water  saturation  S  are  clearly 
shown  in  Fig.  7.  Near  the  entrance,  the  liquid  water  saturation  is 
slightly  higher  than  the  case  with  a  lower  Rch.  But,  as  the  reactant 
gas  moves  downstream,  relative  to  that  with  conventional  uni¬ 
form  flow  field,  the  liquid  water  saturation  with  a  tapered  flow 
field  design  is  lower.  The  liquid  water  saturation  decreases  with 
the  decrease  in  Rc h.  This  implies  that  the  use  of  the  tapered  flow 
field  can  effectively  resolve  the  problem  of  the  water  manage¬ 
ment.  In  addition,  the  liquid  water  saturation  on  the  cathode  side 
is  higher  than  that  on  the  anode  side.  This  is  due  to  the  liquid 
water  generation  occurring  in  the  catalyst  layer  on  the  cathode 
side. 


Fig.  8.  Effects  of  Rc h  on  the  local  axial  velocity  distributions  at:  (a)  V-0.1  V; 
and  (b)  V=0.2V. 


86 


H.C.  Liu  et  al.  /  Journal  of  Power  Sources  158  (2006)  78-87 


In  order  to  explain  the  effectiveness  of  the  tapered  channel  in 
the  water  management,  the  local  velocity  distributions  for  dif¬ 
ferent  7?ch  are  presented  in  Fig.  8.  Here  the  local  velocity,  v(Y), 
is  the  F-distribution  of  the  axial  velocity  at  the  center  of  the  flow 
channel.  It  is  obviously  that  the  velocity  is  higher  in  a  tapered 
channel  with  a  smaller  Rch,  and  the  effect  becomes  strong  with 
the  decrease  in  the  value  of  Rc h.  Therefore,  the  removal  rate  of 
the  liquid  water  is  much  more  effective  for  a  tapered  channel  with 
a  smaller  Rc h.  Moreover,  there  is  a  sharp  increase  of  velocity, 
v,  at  Y>  0.070  since  the  fuel  flow  forced  into  the  porous  layer 
must  be  drained  out  to  channel  at  Rc\l  =  0A.  For  this  reason, 
the  liquid  water  saturation  S  at  the  GDL-catalyst  layer  interface 
decreases  slightly  at  0.070  <  Y<  0.075.  But  at  Y>  0.075  on  the 
GDL-catalyst  layer  interface  near  flow  exit,  the  liquid  water  sat¬ 
uration  S  increases.  This  may  be  due  to  the  fast  accumulation  of 
the  liquid  water  from  the  upstream  and  the  porous  layer,  which 
in  turn,  cause  high  liquid  water  saturation  S. 

In  analysis  of  the  PEM  fuel  cell,  the  pressure  drop  on  the 
cathode  side  is  one  of  the  key  operating  parameters  to  the  cell 
efficiency.  To  examine  this  effect,  Fig.  9  shows  the  results  of 
pressure  drop  on  the  new  flow  channel  design  with  different 
Rch.  It  is  interesting  to  note  that  the  pressure  drops  for  the  cases 
of  tapered  channels  with  Rch  =  0.5  and  0.1,  respectively,  are 
about  2  and  20  times  of  the  uniform-depth  channel  with  Rc h  =  1 . 
Although  there  seems  a  considerable  difference  between  the 


Fig.  9.  Effects  of  Rc h  on  the  pressure  drop  between  gas  fuel  inlet  and  outlet  in 
the  channel  on  cathode  side. 


pressure-drops  at  Rch  =  0.1  and  Rch  =  1,  the  value  of  pressure- 
drop  for  7?ch  =  0.1  is  only  as  small  as  0.018  bar  in  the  typical 
cases  considered  in  Fig.  9.  That  implies  a  small  penalty  with  the 
application  of  the  tapered  fuel  channels. 

5.  Concluding  remarks 

With  an  appropriate  design  of  the  flow  channel  in  the  bipo¬ 
lar  plate,  the  thermal  and  water  management  of  the  PEM  fuel 
cells  can  be  achieved  efficiently.  To  reach  this  end,  a  new 
design  of  tapered  channel  with  non-uniform  depth  has  been  pro¬ 
posed.  The  performance  of  the  tapered  flow  fields  have  been 
investigated  by  using  numerical  predictions  of  the  detailed  gas 
transport  phenomena  and  cell  performance  of  the  PEM  fuel 
cells  with  tapered  flow  field  design.  What  follows  are  the  major 
findings. 

1.  Generally,  the  application  of  tapered  flow  channel  design 
in  the  bipolar  plate  of  PEM  fuel  cells  leads  to  a  favorable 
consequence  in  fuel  transport,  water  management  and  cell 
performance.  The  effects  can  be  enhanced  with  a  decrease  in 
taper  ratio  of  the  fuel  channel. 

2.  The  beneficial  effects  of  the  tapered  flow  channel  become 
remarkable  at  the  conditions  of  low  operating  voltages  or 
high  current  densities.  Therefore,  the  use  of  the  tapered  flow 
channel  will  force  more  reactant  gas  into  the  gas  diffuser 
layer  and  catalyst  layer,  which  in  turn  results  in  a  better  cell 
performance. 

3.  It  is  inevitable  that  a  penalty  of  pressure-loss  accompanies 
the  use  of  the  present  tapered  channels  in  flow  field  design  of 
the  bipolar  plate.  However,  the  present  results  disclose  that 
the  pressure  loss  is  quite  low  in  general  and  the  potential  in 
practical  application  of  this  new  flow  field  design  is  demon¬ 
strated. 

4.  Finally,  the  deviations  between  the  predictions  with  and  with¬ 
out  consideration  of  the  liquid  water  effects  reveal  the  influ¬ 
ences  of  the  liquid  water  on  the  transport  phenomena  and 
performance  of  PEM  fuel  cells.  The  present  results  provide 
evidences  for  significance  of  the  liquid  water  effect  and  imply 
the  requirement  of  inclusion  of  the  effect  in  the  modeling  and 
analysis. 
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